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application to LaFeAsO\_ x F x 



OO 
O 

o 

(N 

as: 
33 

< 
o 

(N 



Yi Zhou 1,2 , Wei-Qiang Chen 1 and Fu-Chun Zhang 1 
Department of Physics, Center of Theoretical and Computational Physics, 
The University of Hong Kong, Hong Kong, China 
Department of Physics, Chinese University of Hong Kong, Hong Kong, China 

(Dated: August 20, 2008) 

We use group theory to classify the superconducting states of systems with two orbitals on a 
tetragonal lattice. The orbital part of the superconducting gap function can be either symmetric 
or anti-symmetric. For the orbital symmetric state, the parity is even for spin singlet and odd for 
spin triplet; for the orbital anti-symmetric state, the parity is odd for spin singlet and even for spin 
triplet. The gap basis functions are obtained with the use of the group chain scheme by taking into 
account the spin-orbit coupling. In the weak pairing limit, the orbital anti-symmetric state is only 
stable for the degenerate orbitals. Possible application to iron-based superconductivity is discussed. 

PACS numbers: 74.20.Rp 
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I. INTRODUCTION 

Symmetry plays an important role in the study of su- 
perconductivity. By using the symmetry of the supercon- 
ducting (SC) gap function, Ginzburg-Landau theory can 
be constructed and electromagnetic response and topo- 
logical excitations can be inspected. In the past decades, 
the symmetry analyses to classify unconventional SC 
states have been focused on single-band superconductors, 
and have shed much light on our understanding of heavy- 
fermion and ruthenate superconductors 

Very recently, a new class of iron-based high tem- 
perature superconductors has been discovered with T c 
as high as above 50K@, i, H i, i, SHU- Ex- 
perimentally, spin density wave (SDW) order has been 
observed in the parent compound LaFeAsO, but van- 
ishes upon fluorine doping where the superconductivity 
appears [13, El- Specific heat measurement as well as 
nuclear magnetic resonance suggested line nodes of the 
SC gap[l2l [l3l [LH. [TH . The transition temperature esti- 
mated based on the electron-phonon coupling is low, and 
unlikely to explain the observed superconductivity [lq] . It 
has been proposed that the superconductivity is of mag- 
netic origin and is unconventional. Local density approx- 
imation (LDA) shows that iron's 3d electrons dominate 
the density of states near the Fermi surfaces in the par- 
ent compound LaFeAsO [HI H d E3, HI El- In their 
calculations, there are three hole-like Fermi surfaces cen- 
tered at the r point and two electron-like Fermi surfaces 
around the M point. By F-doping, the area of the three 
hole-like Fermi surfaces shrinks while the area of the two 
electron-like Fermi surfaces expands. The band struc- 
ture obtained from the LDA may be well modeled by a 
tight-binding model with two or three orbitals (d xz , d yz 



and d xy ) [23, |2j |2J, |2J, |2J, |28|, |29|, |30( . Because of the 
multiple orbitals in the low energy physics, it is natu- 
ral to raise the question how to generalize the symmetry 
consideration from single-band to multi-band cases. 

In this paper, we will generalize the symmetry analyses 
developed for the single band SC state to systems with 



two orbitals. We will use group theory to classify the 
allowed symmetry of the gap functions of the two-orbital 
SC state on a tetragonal lattice by including a spin-orbit 
coupling between the paired electrons. While our focus 
will be on the Fe-based compounds, some of our anal- 
yses may be applied to more general systems with two 
orbitals. 

We arrange this paper as the follows. In Sec. II, we dis- 
cuss the symmetries governing the system and how these 
symmetries affect the Hamiltonian and gap functions. In 
Sec. Ill, we consider the possible two-orbital SC states 
on a tetragonal lattice. Section IV is devoted to sum- 
mary and discussions. We also supply some appendixes 
for details. In Appendix \Kl we show how the symme- 
tries give rise to the requirements to the non-interacting 
Hamiltonian. In Appendix [Bl we specify the point group 
D^h of lattice according to space group PA/nmm. In 
Appendix [C] we discuss how the gap functions transfer 
under symmetry operations. In Appendix IDl we discuss 
the energy gap functions in the degenerate bands. 



II. SYMMETRY OF GAP FUNCTION A (k) 

We consider a tetragonal lattice, appropriate for doped 
LaFeAsO. Since our primary interest is in the SC state, 
we will not consider the translational symmetry broken 
state such as the spin density wave state observed in the 
parent compound of LaFeAsO. The system is invariant 
under both time reversal and space inversion. The in- 
version symmetry suggests that the SC pairing is either 
even or odd in parity. We shall assume in this paper that 
the time reversal symmetry remains unbroken. 

We consider a system described by Hamiltonian 



H = H + Hp, 



(1) 



where Hq is non-interacting part, and H pa i r is a pairing 
Hamiltonian, and H so is the spin-orbit coupling of the 
Cooper pairs. We shall consider the SC state preserves 
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all the symmetries in Ho except the U(l) symmetry in 
electric charge and the spin rotational symmetry due to a 
weak H so . We assume H to be given by a tight-binding 
Hamiltonian 

Hq = y ] cj CQiS ^l C QiQ 2 C kQ2Sj (2) 

where a = 1,2 are the orbital indices, which correspond 
to the two orbitals 3d xz and 3d yz in Fe, s =|, | are the 
spin indices. Note that for LaFeAsO, the actual crys- 
tal structure has two Fe-atoms in a unit cell due to the 
As atomic positions, which are allocated above and be- 
low the Fe-plane alternatively. For convenience, here we 
use the extended Brillouine zone, and the summation k 
is in the extended zone. Hq is invariant under symme- 
try transformation. This requires certain symmetries on 
£kaic(2! which we will discuss in detail in Appendix [A"l 
We note that a more appropriate model should also in- 
clude the da;j,-orbital[25], but we shall leave the symmetry 
analyses of the three orbitals for future study, and con- 
sider a simplified version of the two orbital case in this 
paper. 




(a) (b) 

FIG. 1: (color online) Lattice structure of LaFeAsO. It is a 
tetragonal lattice with two Fe atoms per unit cell. The lattice 
constants are a — b ~ 4.034 and c ~ 8.74A[li"l]. where a is 
the distance between two next nearest neighbor Fe atoms, 
(a) Origin choice 1 of space group PA/nmm, at 4m2 and at 
(—a/4, a/4, 0) from center (2/m). It can be chosen either at 
an Fe or at an O atom; (b) Origin choice 2 of space group 
PA/nmm, at center (2/m) and at (a/4, —a/4,0) from 4m2. It 
can be chosen either at the midpoint of two nearest neighbor 
Fe atoms or at the midpoint of two nearest neighbor O atoms. 
Here 2/m denotes the two fold rotation C2 and reflection m 
(see Appendix iBl for details). The origin choice 1 and 2 are 
different from each other by a shift of f— a/4, a/4. 0l[3lT]. 



The gap function of the two-orbital SC state can be 



generally written as 

^si 1 ^ 2 (k) = — V£V?s"s" i (k) k ) ( c k'a 3 s 3 C-k'a4S4), 

k'a3CK4 

S3S4 

(3) 

where V^"|"g "^ Q4 (k, k') is the effective attractive inter- 
action. Hereafter we will use the matrix notation A (k) 
for the gap function. 

To classify the symmetry of the SC gap function for 
multiple orbitals, we recall that in the single orbital case, 
the spin-orbit coupling of the Cooper pair plays an im- 
portant role to the non s-wave superconductors, and the 
symmetry of the gap function is determined by the crys- 
tal point group of the lattice and the spin part of the gap 
function. In the two-orbital system, the orbital degree 
of freedom is usually coupled to the crystal momentum, 
hence to the spin via the spin-orbit coupling. Therefore, 
the spin, spatial, and the orbital parts are generally all 
related in the gap function. 

Let us first discuss the crystal symmetry. The crystal 
structure of LaFeAsO is shown in Fig. 1. The tetrago- 
nal crystal symmetry is characterized by the point group 
D±h. The tetragonal point group may be specified ac- 
cording to the space group P4/nmm of the compound, 
and the details will be discussed in Appendix [Bl There 
are five irreducible representations of D4 group, denoted 
by r, including 4 one-dimensional representations (A\, 
A2, Bi and B2) and 1 two-dimensional representation 
(-E^dH. The tetragonal lattice symmetry requires H to 
be a "scalar" or A\ representation of D4. In the absence 
of spin-orbit coupling, spin is rotational invariant and we 
have both the point group symmetry and the spin rota- 
tional symmetry. [32| 

We now discuss the orbital degrees of freedom in con- 
nection with the crystal symmetry. The two orbitals d xz 
and d yz transform as E representation of -D4. In general 
the orbital indices d xz and d yz are not good quantum 
numbers because of the mixed term of the two orbitals 
in Ho, and the two energy bands are not degenerate. In 
that case it is necessary to include the coupling of the 
orbital to spatial and spin degrees of freedom. 

Without loss of generality, the gap function can be 
written as a linear combination of the direct products of 
the orbital part Vl and the spin part A spm in a given 
representation T of the point group D4, 

A(T;k) = V(r,m) (T,m\T LS ,m LS ;Tn,m n ) 

m,r LS ,rsi 

xA s P m (r LS ,m iS ;k)<^(r n ,mo), (4) 

where both A spm and are 2x2 matrices, A^™ dic- 
tates the pciiriiig in spin sps>C6 cmd i^o^ 0.2 

dictates the 

pairing in orbital space, T^s and Tq are irreducible rep- 
resentations of D4 in spin and orbital spaces, respectively, 
ra,mis,ra[j are bases of representations r, Tls,Tq, re- 
spectively, (r, m\TLs, ttilsI Ffi, ma) is the Clebsch- 
Gordan (CG) coefficient. Note that the k-dependence 
is contained in A spm , but not in fl. Here ry (r, m) is the 
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coefficient of the basis m of the representation T. The 
anti-symmetric statistics of two electrons requires 

A T (-k)=-A(k). (5) 

Below we will first discuss A s ? m (r L5 ;k) and fi (r n ) 
separately, and then combine the two to form an irre- 
ducible representation T of D4. We follow Sigrist and 
UedaQ and write A s *> m (k) in terms of the basis func- 
tions ip (r, to; k) for the spin singlet S — and d (T, to; k) 
for the spin triplet 5=1, 

A spm (r, m;k) = i [a *jj (T, m; k) + a ■ d (T, m; k)] a 2l 

(6) 

Here ip (k) is a scalar and d (k) is a vector under the 
transformation of spin rotation. For this reason, it 
is more convenient to use ip (k) and d (k) instead of 
A spm (k) to classify the pairing states. 

Due to the fermionic anti-symmetric nature, the gap 
function must be anti-symmetric under the two par- 
ticle interchange, or under a combined operations of 
space inversion, interchange of the spin indices and in- 
terchange of the orbital indices of the two particles. 
Let Pi,2 be the two particle interchange operator, and 
Pspace, Pspin, Porbital be the interchange operator acting 
on the space, spin, and orbital, respectively, then the 
fcrmion statistics requires 

P\,2 — P s P ac&P spinP orbital L (7) 

Since the system is of inversion symmetry, the pairing 
states must have either even parity P spa ce = +1 or odd 
parity P spa ce = — 1- Furthermore, the total spin S of the 
Cooper pair is a good number, and this is so even in the 
presence of H so , which breaks spin rotational symmetry 
but keeps inversion symmetry, so that it does not mix 
the S = 1 with S = states. Therefore, under the two 
particle interchange, the spin part of the gap function 
must be either symmetric: (P sp i n = +1, with S = 1) 
or anti-symmetric (P sp i n — — 1 with S = 0), represented 
by the vector d or the scalar if) in Eq. ([6]), respectively. 
Because of the inversion and spin symmetries, we have 

Porbital — il. 

The orbital part of the pairing matrix is spanned 
in the vector space of (d xz ,d yz ), which is an irreducible 
representation E of the point group D±. Thus f2 be- 
longs to an irreducible representation given by E ® E = 
A\ © A2 © B\ © B2, which are all one-dimensional, hence 
simplifies the classification of the pairing states. Accord- 
ing to the CG coefficients of D4 group, up to a global 
factor, = Co in representation A±, = 03 in Bi, 
and f2 = o\ in B2, which are all orbital symmetric: 
Porbital = +1. = a% in A2 representation, which 
is orbital anti-symmetric: Porbital = — 1. In brief, A\ 
and B\ of f2 are representations for intra-orbital pair- 
ing, B2 is for symmetric inter-orbital pairing and A2 is 
for anti-symmetric inter-orbital pairing. For convenience, 
we choose fi to be Hermitian, so that ip (k) and d (k) 
will be real. (Wan and Wang[36] pointed out that Pauli 



matrices transfer as four one-dimensional irreducible rep- 
resentations.) 

The crystal point group of the lattice will dictate the 
allowed symmetry in k space. The transformation of ip 
and d under symmetry operations can be found in Ap- 
pendix [UJ In Sec. Ill, we will study the basis functions 
"0(r,TO;k) and d(r,m;k), and combine them with the 
orbital part Q to obtain the irreducible representations 
of group D4. 

III. POSSIBLE TWO-ORBITAL SC STATES ON 
A TETRAGONAL LATTICE 

We will use the group chain scheme to study the repre- 
sentation and the basis function of ip and d by assuming 
a spin-orbit coupling. In the group chain scheme, we be- 
gin with a rotational invariant system in both spin and 
spatial spaces. The representation of its symmetry group 
G can be decoupled into a spatial part Dr^s and a spin 
part D(s)i with L as the relative angular momentum of 
the Cooper pair, 

D {G )=D(L)®D {S) , (8) 

In the presence of the spin-orbit coupling, D(l) an d D(S) 
are no longer the irreducible representation of the rota- 
tional group, but the total angular momentum J = L + S 
is, and D( j) is the corresponding irreducible representa- 
tion of the rotational group. 

We now turn on a crystal field with tetragonal lattice 
symmetry group D4, so that the rotational group SO(3) 
is reduced to D4, and D(l) ® D(S) is reduced to a direct 
product of irreducible representations Tj^s of group -D4, 

D (L) ®D (S) (9) 

Including the coupling to the orbital part f2, the repre- 
sentation D(y ls -j <8) D(r,n) is decomposed into irreducible 
representations , 

D(r LS )®D (rn) =@D (r) . (10) 
r 

D(r n ) ls one-dimensional, thus these representations have 
a very simple form. 

Let us consider the even parity case. From Eq. (7), 
the SC gap function can be either orbital symmetric 
Porbital = +1, spin singlet or orbital anti-symmetric 
Porbital = —1, spin triplet. We list the SC gap basis 
functions for spin singlet and spin triplet according to 
the irreducible representations T in Tables I and II re- 
spectively. The listed even pairing states include s-wave 
(extended s-wave), d-wave and g-wetve. Here 0, 0, 2, 2, 
and 1 are natural notation for the five irreducible repre- 
sentations of Dih\ A\, A2, E>i, B2 and E are Schonflics 
notation; Ti_5 are Koster notation. According to Eq. 
©, the gap function of the SC state is a linear combi- 
nation of the basis functions in one irreducible represen- 
tation r, and the basis functions belonging to different 
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representations in T, e.g. A\ g and Big, will not mix with 
each other. 

We are particularly interested in 2D or quasi-2D lim- 
iting cases, relevant to Fe-based SC compounds, where 
the gap function is fc z -independent, and the Fermi sur- 
face is cylinder-like. However, for completeness we also 
list in the Tables those three-dimensional basic functions 
marked with 3D. 

In the last column of each table, we list the allowed 
energy zeroes in the quasiparticle dispersion determined 
by the gap functions for the special case that the two 
energy bands are completely degenerate. The detailed 
calculations for the quasiparticle energies in the degener- 
ate cases are given in Appendix [Pi We will discuss the 
quasiparticle properties for the non-degenerate cases in 
the discussion section below. 

Similarly, for the odd parity pairing P spa ce = — T we 
can have either orbital anti-symmetric P or utai = — 1, spin 
singlet, or orbital symmetric P or utai = 1, spin triplet, 
which are listed in Tables III and IV respectively. For the 
spin triplet, we list p-wave, /-wave and /i-wave states. 



IV. SUMMARY AND DISCUSSIONS 

In summary, we have studied the pairing symmetry 
of the two orbital superconducting states on a tetrago- 
nal lattice. Based on the symmetry consideration, we 
have classified symmetry allowed pairing states with the 
space inversion, spin, orbital, and the lattice symmetries 
by including a spin-orbit coupling. In addition to the 
even parity for the spin singlet and odd parity for the 
spin triplet pairings, familiar in the single band supercon- 
ducting gap functions, which corresponds to orbital sym- 
metric pairing in the two orbital systems, there are also 
even parity for spin triplet and odd parity for the spin 
singlet pairings, corresponding to orbital anti-symmetric 
pairing. The symmetry allowed gap basis functions are 
listed in Tables I- IV in the text. In the orbital sym- 
metric states, the gap basis functions within the same 
representation of the point group but with different or- 
bital representations are allowed to combine to form a 
gap function. 

Below we shall discuss some limiting cases. First, we 
consider the weak pairing coupling limit. In this case, 
we can diagonalize Hq firstly to obtain the two energy 
bands. H pa i r in Eq. JT]) is to induce a pairing of elec- 
trons near the Fermi surfaces within a very small energy 
window. If the two energy bands are not degenerate, then 
the two Fermi surfaces do not coincide with each other, 
and the pairing will only occur between electrons in the 
same band, since the energy mis-match of the two elec- 
trons with opposite momentum in the two bands will not 
lead to the SC instability in the weak coupling limit. The 
issue is then reduced to the two decoupled single band 
problem. Because the intra-band pairing is between sym- 
metric orbitals, all the states with orbital anti-symmetric 
pairings such as those listed in Tables II and III will not 



be realized. There is a one-to-one correspondence be- 
tween the present work and the single band analysis 
In terms of the orbital picture, the intra-band pairing 
gap function is described by a linear combination of the 
orbital representations a®, 01,03 in each representation 

of r. 

The strong pairing coupling case is more complicated, 
and possibly more interesting. The symmetry analyses 
we outlined in this paper may serve as a starting point. 
The pairing interaction may overcome the energy mis- 
match of the paired inter-band electrons to lead to the 
superconductivity. In a recent exact diagonalization cal- 
culation for a two orbital Hubbard model on a small size 
system, Daghofer et al. have found an inter-orbital pair- 
ing with spin triplet and even parity with the gap func- 
tion to be cos k x + cos k y [H[ . Their pairing state corre- 
sponds to Eg representation in Table II, and provides a 
concrete example of the orbital anti-symmetric pairing 
state. Generally we may argue that the gap structure 
will be gapless with Fermi pockets for 2D systems un- 
less the pairing coupling is strong enough to overcome 
all the mis-matched paired electrons in the momentum 
space. An example was given by Dai et al.[23| and also 
discussed by Wan and Wang[3y|. This seems to essen- 
tially rule out any possibility for line nodes in the orbital 
anti-symmetric pairing state in the strong pairing cou- 
pling limit. A nodal in quasi-particle energy requires 
the gap function to vanish. As a result, the pairing 
strength near this nodal will not be strong enough to 
overcome the energy mis-match of the inter-band paired 
electrons. Therefore, a nodal in quasi-particle energy im- 
plies a Fermi pocket in this case. 

Another interesting limit is the two orbitals are com- 
pletely degenerate: £kai,a 2 = Ck^ai,a 2 - The system has 
an orbital SU(2) symmetry. In this case, our analyses are 
most relevant, and all the classified states listed in Tables 
I-IV could be stable even in the weak pairing interaction. 
Because of the orientational dependence of the orbitals 
in crystal, such degeneracies may not be easy to realize. 
A possible realization is on the materials with two-fold 
pseudospin symmetry or two-valley degeneracy such as 
in graphene. While the point group will depend on the 
precise crystal symmetry concerned, but some general 
features discussed in this paper may be applied to those 
systems. 

We now discuss the band structure in the extended 
zone and the reduced zone. Because of the positions of As 
atoms, the translational lattice symmetry is reduced and 
the Brillouine zone is halved. In general, such a transla- 
tional symmetry reduction may lead to hopping matrix 
between momentum k and k + Q in the extended zone, 
with Q = (7r, 7r)/a' and a' — a/y/2 is the lattice con- 
stant of reduced unit cell. However, for the two orbitals 
d xz and d yz , the point group symmetry prohibits the hy- 
bridization between states at k and k + Q, if we only 
consider intra-layer hopping. The tight-binding Hamil- 
tonian adopted by both Raghu et al.[27| and Lee and 
Wen[25[ explicitly illustrate the vanishing of the mixing 
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TABLE I: Superconducting gap basis functions tp(k) on tetragonal lattice for even parity, orbital symmetric and spin singlet 
pairing states. F: representation of D4. The listed notations are natural, or Schonflies and Koster (in parentheses). Q: orbital 
representation, 00 is the identity matrix, and 01,2,3 are Pauli matrices. Listed gaps properties are for the two completely 
degenerate orbitals. fc z -dependent basis functions are marked with (3D), listed for completeness. 



r 


basis ?/>(k) 


Q 


gap 


(Ax g , r+) 


I5 kx fey'-) k z (^-D) 


00 


line nodal, or full gap 


k x " ky 


03 


k x ky 


01 


(A 2g , r+) 


kxky(^k x ky^ 


00 


line,full 


k x ky 


03 


,2 7,2 

™x ™y 


01 


2 (Bi a , r+) 


kl-kl 


00 


line,full 


^-■>k x J rk y ]k z (3D) 


03 


k x ky{k x ^1/) 


01 


2 (B2 3 , r+) 


k x ky 


00 


line,full 


kxky(k x 


03 


I5 k x + k y ] k z (3D) 


01 


1 (E g , r+) 


{kxkzikykz} (3D) 


00, 03, 01 





TABLE II: Superconducting gap basis functions d(k) on tetragonal lattice for even parity, orbital anti-symmetric and spin 
triplet pairing states. Notations are the same as in Table I. 



r 


basis d(k) 


Q, 


gap 


(Alg, r+) 


2, (k x + k y )z, {k x + k y )z, k x k y z 


02 


line, full 


6 [A ag , r+) 


k z (k x y - k y x) (3D) 


02 




2 (Big, r+) 


{hi - k'y)z; k z (k x x - k y y) (3D) 


02 


line 


2 (B 2g , r+) 


k x kyZ\ k z {k x x ~f~ kyij^ (3D) 


02 


line 


1 (E g , r+) 


{%■> y)i {k x %i k x y), [k y x, k y y), {k x k y x, kxk y y)' 7 
(k z Xi k z y^ 7 {k x k z z , kyk z z} (3D) 


02 


line, full 



term. Therefore, we may discuss the SC symmetry using 
the extended zone and using Hq given in Eq. In the 
extended zone, there is only one Fermi point for each k, 
hence the bands are not degenerate. In the weak pair- 
ing coupling limit, all the orbital anti-symmetric pairing 
states will be irrelevant, and the weak coupling theory 
will naturally lead to the orbital symmetric states. 

Near the completion of the present work, we learned 
of the similar work by Wan and Wang. [36| , who consid- 
ered SC symmetry for two-orbital pairing Hamiltonian. 
Our results are similar to theirs, with the difference that 
we have included a spin-orbit coupling term in our group 
theory analysis, while this term was not explicitly in- 
cluded in Ref . [3q . As a result, our classification for the 
spin triplet states is not the same as theirs. Such differ- 
ence may be amplified when we discuss some behaviors 
related to spin degrees of freedom. We also note that 
similar group theory analysis were carried out for the 
two band pairing Hamiltonian by Wang et al.[34j. Since 
they adopted the bands instead of the orbitals, a direct 
comparison is not apparent. We become aware of another 
related work[33| after completing the present work too. 
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APPENDIX A: THE SYMMETRY OF £ kQ1Q2 IN 
EQUATION ([2]) 

In this appendix, we will discuss the symmetry re- 
quirement of £kaia 2 - The non-interacting Hamiltonian 
given by Eq.([2]) should keep invariant under any sym- 
metry transformation of point-group D4, hence Hq be- 
longs to the representation A±. This symmetry Re- 
quirement will affect the choice of £kaia 2 - F° r con- 
venience, we use the 2x2 matrix form £k in orbital 
space, thus £k can be rewritten in terms of Pauli ma- 
trices, |k = Ck CT o + £ k ci + £k CT 2 + £ k CT 3- Similarly to the 
case of ft, (pl s a (/)us, <t> ks cri<j>k s , (frts^ks and ^ ks 
transform as A±, B2, A2 and B\ respectively, where 
4>ks = (ckr, Ck2) T - Using the CG coefficients of D4 — C4 
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TABLE III: Superconducting gap basis functions ?/>(k) on tetragonal lattice for odd parity, orbital anti-symmetric and spin 
singlet pairing state. Notations are the same as in Table I. 



r 


basis y)(k.) 




gap 


(Aiu, IT) 


k z (3D) 


£>"2 




o (A 2u , it) 


k z (ki - 6k 2 x k y + fcj) (3D) 


C"2 




2 r-) 


fc«(*2 - *g) (3D) 


0"2 




2 (B 2 u, rj) 


k x kyk z (3D) 


0"2 




i (Eu, r-) 


(/l X , ky ) 


0"2 


line 



TABLE IV: Superconducting gap basis functions d(k) on tetragonal lattice for odd parity, orbital symmetric and spin triplet 
pairing states. 



r 


basis d(k) 


ft 


gap 


(Am, IT) 


fc^a; + k y y; k z z (3D) 


O"0 


line,full 


fc>xx fcyy 


CT3 


kyX 


ui 


6 (A 2u , IT) 




O"0 


line,full 










2 (Biu, IT) 




Co 


line,full 


k x x + fc a y; fe z i (3D) 




kyX kxy 


ui 


2 (s 2u , rj) 


kyX -{" fcx^ 


O"0 


line,full 


kyX k x y 


0"3 


k x x + k y y; k z z (3D) 


o"i 


i (Eu, r 5 ") 


(k x z,k y z) ;(k z x,k z y) (3D) 


a a , cr 3 , a i 


line 



group chain{33j, we find that £ k > and transform 
as i?2, ^2 and i?i respectively. Some examples of 
£ k ' ' ' are shown in the following, 

£k = 1, cos k x + cos k y , cos k x cos k y , 
= sin fca; sin fcj, , 

= sin fca; sin fc y (cos k x — cos fcy) , 

£k = COS /Ca; — cos k y . 

APPENDIX B: THE POINT GROUP D ih 

Here we would like to specify the tetragonal 
point group according to the LaFeAsO space group 
P4/nmm.[3l| In real space, the point group is nei- 
ther usual D^h = D4 ® ah nor usual -D4 generated by 
{CT, C2 y }, where <ih is the reflection refer to xy plane, 
Ci z is the four-fold rotation around the z axis, and Gi v is 
the two-fold rotation around the y axis, where (x, y, z) is 
specified in Fig. [51 However, it contains two subgroups, 
which refer to two different origin choices of the lattice. 
One is a subgroup of generated by {C^tx^, Ci y Oh\, 
which is also a D<± group (or to be precise, Did, an iso- 
morphic group to Di) with origin choice 1 as shown in 



Figs, [lja) and[UJa). The other subgroup is the direct 
product of inversion symmetry group / and cyclic group 
C2xy with origin choice 2 as shown in Figs. [Tfb) and 
[2f b) . The transformation of (x, y, z) under these symme- 
try operations can be found in Tables [V] and IVII Hence 
in k-space, it is still a tetragonal point group D^h- 



O O OFe O O O 




• • • • 

o o o o o o 
® <b) 

FIG. 2: (color online) FeAs layer and specification of (x,y, z). 
xy plane consisting of Fe atoms is shown in (a) and (b) with 
different origin choices, z axis is perpendicular to Fe-plane. 
"+" denotes an As atom above the Fe-plane, while "—"denotes 
an As atom below the Fe-plane. (a) Origin choice 1. (b) 
Origin choice 2. (also see Fig[TJ 

There are five irreducible representations of group, 
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TABLE V: Eight symmetry operations of D2d (an isomorphic 
group to Di) will generate eight general positions. Where 
"general" is defined as the following: a set of symmetrical 
equivalent points is said to be in "general position" if each 
of its points is left invariant only by the identity operation 
but by no other symmetry operation of the space group. The 
origin choice is 1. 



group element 


general position 


E 


(x,y,z) 




(y, -x, -z) 


{CAzOhf 


(-x,-y,z) 


{CAzO h ) 3 


{-y,x, -z) 




(-x,y,z) 


C2yUhC4zO'h 


{-y,-x,-z) 


C'iyOh (CAzOhj 2 


(x,-y,z) 


ClyOh {CAzOhj A 


(y,x, -z) 



TABLE VI: Four symmetry operations and corresponding 
general positions of group / x C-2xy- The origin choice is 2. 



group element 


general position 


E 


(x,y,z) 


d 


{-x, -y, -z) 


Czxy 


{y,x,z) 


GiC^Zxy 


{-y, -x, -z) 



four of them, A%, A2, B\ and B2 are one-dimensional 
representations, and one of them, E is a two-dimensional 
representation. All these five representations are rep- 
resentations of group Dih too. However, there are 2 
two-dimensional irreducible representations E' and E" 
of Dih group, neither of them is the representation of 
group D4. Naively, E' and E" can be viewed as subrep- 
resentations of two irreducible representations of group 
SU(2), J = 1/2 and J = 3/2, respectively. Since the 
representations E 1 and E" can not result in quadratic 
terms in Hamiltonian or Ginzburg-Landau free energy, 
we will not discuss them in this paper. 



APPENDIX C: TRANSFORMATION OF GAP 
FUNCTIONS 

It is not A (k) but ip (k) and d (k) transform as rep- 
resentations of symmetry group. In this appendix, we 
list the transformations of ip (k) and d (k) under vari- 
ous symmetry operations. Firstly, under a point-group 
transformation g, ip (k) and d (k) transform as 

3 V(k) = ^(-D(- G) (p)k), 

gd(k) = D+jGrtd^Gtfk), (CI) 

where D ^ (g) is the representation in three-dimensional 
space with positive (spin-space) or negative (k-space) 



respectively. Secondly, time-reversal transformations of 
ip (k) and d (k) take the forms, 

Kip (k) = ip* (-k) , Kd (k) = -d* (-k) . (C2) 

The anti-symmetric nature of Fcrmion systems, see Eq. 
©, will lead to 

V(-k) = V(k),d(-k) = -d(k), (C3) 

for symmetric and 

^(-k) = -^(k),d(-k)=d(k), (C4) 

for anti-symmetric Q. Hence, combining the above and 
the Hermitian choice of f2's, the time-reversal invariance 
conditions for if> (k) and d (k) become 

r(k)=V(k),d*(k) = d(k), (C5) 

since under time-reversal transformation, Q transforms 
as 

Kfl = ft*. (C6) 

APPENDIX D: ENERGY GAP FUNCTIONS IN 
THE DEGENERATE BANDS 

The energy gap of the superconducting states indeed 
depends on the details of interaction, especially, depends 
on the ratio of St /A, where 5t is the energy scale of the 
splitting of two bands and A is the energy scale of pair- 
ing potential. In the the "strong pairing coupling" limit 
St <C A, we expect the energy gap is close to dt = 
case, say, two bands are degenerate. A small perturba- 
tion proportional to St/X will not change the energy gap 
very much, e.g. close the full gap or change from full 
gap to line nodal gap. In the weak pairing coupling limit 
A <C St, the situation may be very different from strong 
coupling limit, which is discussed in Ref . (36| . So that we 
will focus on the strong coupling limit and assume two 
degenerate bands in the following. 

Due to two degenerate bands, the effective mean field 
Hamiltonian in k-space can be written as an 8 x 8 matrix, 

fr _( Ckoo ® cr A(k) \ , . 

Hk -{ At(k) -^o^oj' (D1) 

with the basis Ck = 

(ckti, c k |2, Ckii, Ckj.2, cL kn , cL kT2 , cL kil , cl ki2 ) T . The 
indices in the 4x4 matrices A(k) and £kOo ® a o are 
arranged as the following, by direct products the former 
two indices denote spin space and the later two denote 
two orbitals. It is easy to know the energy dispersion, 

Ei,, = ±v/£ + A£ M , (D2) 

where A k/J is one of the eigenvalues of the matrix 
A(k)At(k). For degenerate bands, the minimum of |A k/ J 
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is the energy gap. For simplicity, we will focus on the k z - 
independent pairing with a cylinder-like Fermi surface 
which is the case of doped LaFeAsO most likely. 

At first, we will consider the even parity, orbital anti- 
symmetric, spin triplet pairing states in Table II. Or- 
bital anti-symmetric states have only one component a 2 
in the part. Gap function is of the form, A(k) = 
i [a ■ d(k)]<7 2 <g> 02, thus Aj^ = |d| 2 ± |d x d*|. For the 
time-reversal invariant state, d = d*, the gapless condi- 
tion follows as |d| 2 = 0. For B\ g and B 2g states, they 
are d-wave states and have line nodal gap at Fermi sur- 
faces. A\ g states can be of either s-wave or extended s- 



wave. The s-wave state is of full gap while the extended 
s-wave state possibly has line nodal gap at Fermi surface, 
e.g., the state d (k) = cosk x cos k y z. The E g representa- 
tion involves s-wave, extended s-wave and d-wave states. 
The s-wave state is fully gapful, the d-wave state has line 
nodal gap, the extended s-wave state can be either fully 
gapful or of line nodal gap. 

Then we consider the odd parity, orbital symmetric, 
spin triplet pairing states in Table IV. Orbital symmetric 
states have three components cro,i,3 m the f2 part. Gap 
function can be written as 



A(k) = i [a ■ d (k)] <r 2 ® cr + i W ' d i ( k )] °2 ® v\ + i [<J ■ d 3 (k)] a 2 ® a 3 , (D3) 

thus 

A(k)At(k) = [(|d | 2 + |d 1 | 2 + |d 3 | 2 ) ( 7o+2(doxd;+d 1 xdt+d 3 xd*)-a]® ( 7o 
+ [(d • d{ +di • do) CTo +«(d xdi +djxdo) ■ a] ® ai 
+ [(d • d 3 + d 3 • do) cj o + i (d xd* 3 + d 3 xd^) • a] ® a 3 

+ [(dixd* - d 3 xdl) ■ g - i (di • d* - d 3 • dj) <T ] <g> a 2 . (D4) 
For a time-reversal invariant state, d* (k) = dj (k), i = 0, 1, 3, so that the above can be simplified as 

A(k)A f (k) = (d 2 , + d 2 + d§) CT ® cr + 2 (d • di) a a ® n + 2 (d • d 3 ) a ® cr 3 + 2 (dixd 3 ) • a ® a 2 . (D5) 
We obtain from the above 



A 2 ^ = (d 2 + dj + df) ± 2yJ (d • di) 2 + (do • d 3 ) 2 + (d lX d 3 ) 2 . (D6) 

Gapless condition reads 

d 2 + d 2 + d 2 - 2\j (do • di) 2 + (d • d 3 ) 2 + (d lX d 3 ) 2 . (D7) 



Careful analysis shows node can appear only when at 
least one of |do|, |di| and |d 3 | vanish. So that E u states 
in Table IV are of line nodal gap. The other four repre- 
sentations A\ u , A 2u , B iu and B 2u can be of either line 
nodal or full gap. For example, for an A lu states in Ta- 
ble IV which consists of two components in the f2 part, 
d (k) = sin k x x + sin k y y, d 3 (k) = sin k x x — sin k y y and 
di (k) = 0, nodal lines will appear at sink x — and 
sin k y — 0. Moreover, any A\ u state in Table IV which 



consists of only one component in the f2 part is of full 
gap. 

Similar consideration will lead to the results for spin 
singlet states shown in Tables I and III. Of course, when 
the ratio dt/X becomes large, the situation may change. 
This change strongly depends on the details of both pair- 
ing states and the Hamiltonian. For example, for an s- 
wave with d (k) = z, and fl = a 2 , Fermi pockets may 
appear when St and A are comparable. 
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